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^ ! Abstract 

f~^ . We derive new constraints on the products of explicitly i?-parity violating couplings 

^> I A and A' in MSSM from searches for /i-e conversion in nuclei. We concentrate on the 

rS^ ' loop induced photonic coherent conversion mode. For the combinations |AA| which 

^~^- in fi-e conversion can be probed only at loop level our constraints are in many 

Qh! cases more stringent than the previous ones due to the enhancement of the process 

(— I I by large ln(m^/m^). For the combinations of |A'A'| the tree-level //-e conversion 

IJ ' constraints are usually more restrictive than the loop ones except for two cases which 

involve the third generation. With the expected improvements in the experimental 
sensitivity, the /U-e conversion will become the most stringent test for all the involved 



X 



C^ ' combinations of couplings. 



1 Introduction 



The minimal supersymmetric standard model (MSSM) extended by explicit 
breaking of -R-parity, R = (— 1)'^(-^~^)+2'5' [xj^ j^as recently received a lot of 
attention. In MSSM, the conservation of i?-parity is put in by hand. Therefore, 
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it is theoretically not well motivated. On the other hand, models exist in 
which the violation of -R-parity is a necessity, e.g., the supersymmetric left- 
right model [2]. Furthermore, the observed excess of high Q^ events in HERA, 
if interpreted in the framework of supersymmetry, seems to indicate non-zero 
-R-parity violating couplings of squarks [3]. 

Limits on i?-parity violating couplings have been derived from a large variety 
of observables. Since the presence of -R-parity violation leads automatically to 
the presence of lepton and/or baryon number violating processes, searches for 
them are especially suitable for constraining these couplings. In this Letter we 
derive new upper limits on the products of the couplings Ajj^ and A^^-^ from 
the non-observation of /i-e conversion in nuclei occuring at one-loop level. We 
compare our results with the previously obtained constraints including the 
ones derived from the tree-level /x-e conversion. The latter process provides 
particularly restrictive bounds on some products of the couplings [4] , but not 
for all of them. The loop induced conversion although suppressed by the loop 
factors is partly sensitive to different couplings than the tree-level process. 
Moreover, it was pointed out in Ref. [5] that in a wide class of models with 
effective interactions of four charged fermions, thus including the MSSM ex- 
tended by -R-parity violation, the /i-e conversion at one-loop is enhanced by 
large logarithms. Therefore, as will be seen, the /i-e conversion at loop level is 
more sensitive to certain combinations of couplings than the tree-level process. 

The motivation for the present work is twofold. Firstly, our study enables us to 
put more restrictive bounds on some combinations of the couplings than do the 
previous works. Secondly, the expected improvements in the sensitivity of the 
/x-e conversion experiments running presently at PSI will make /i-e conversion 
the main probe of muon flavour conservation for most of the extensions of the 
standard model allowing for an order of magnitude more stringent tests of 
-R-parity violation already in forthcoming months. Indeed, while the present 
bound on the conversion branching ratio in 22^^ is -R^* ~ 4.3 ■ 10~^^ [6], the 
SINDRUM II experiment taking currently data on gold, ^^ Au^ should reach 
-R^" ^ 5 ■ 10^^^ and next year on titanium -R^* ^ 3 ■ 10~^^ [7]. Furthermore, 
very recently a proposal for an experiment called MECO has been submitted 
to BNL [8] which is planned to achieve the sensitivity better than 10^^^. 
The experimental prospects for searches of muon number violation have been 
recently reviewed by A. Czarnecki [9]. 



2 Previous constraints 



To compare our analyses with the previous ones, we first collect and list the 
updated upper bounds on all products of -R-parity violating couplings which 
are testable in /x-e conversion. Within the MSSM particle content the gauge 





previous bounds 


//-e at loop level/ \/B 


m? = 
100 GeV 


100 GeV 


1 TeV 


-^121 -^122 
^131 -^132 
A231 A232 
A231 Ai31 
A232 Ai32 
A233 A133 


6.6 • 10-7 [30] 
6.6 • 10-7 [30] 
5.7-10-5 [29] 
6.6 • 10-7 [30] 
1.1 • 10-4 [29] 
1.1 • 10-4 [29] 


4.2 • 10-6 
5.3-10-6 
5.3-10-6 
8.4 - 10-6 
8.4 - 10-6 
1.7- 10-5 


3.2 - 10-4 
3.9 - 10-4 
3.9 - 10-4 
6.4 - 10-4 
6.4 - 10-4 
1.0- 10-3 



Table 1 

Upper limits on the products |AA| testable in /U-e conversion for two different scalar 
masses m? =100 GeV and nif =1 TeV. The previous bounds scale quadratically 
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with the sfermion mass. The scaling factor B is defined in the text and currently 
B = 1. 

invariance and supersymmetry allow for the following i?-parity violating su- 
perpotential [10] 

Wf, = \,uhL,El + K.MPl + K.kU^D'pi - fi,UH, , (1) 



where Xijk = —Xjik and X'^jj^ = —X'lkj- The A, A' and /x terms violate the 
lepton number, whereas the A" terms violate the baryon number by one unit. 
The last bilinear term in Eq. (1) gives rise to interesting physics which has 
been studied elsewhere [11]. We will not consider these effects in the present 
paper and we take as our working model the MSSM extended by the trilinear 
i?-parity violating terms in Eq. (1). 

Simultaneous presence of L violating (Xijk, X'ijk) and B violating {X'-jiJ cou- 
plings would lead to too fast proton decay unless |A' ■ A"| < IQ-^ [12] for 
squark masses below 1 TeV. Several examples exist of models, in which there 
are huge differences between the strengths of lepton and baryon number vi- 
olating couplings. In models, which break i?-parity spontaneously, only the 
lepton number is violated. There are also examples of GUT models, in which 
quarks and leptons are treated differently [13,14]. In the following we assume 
that only lepton number violating couplings are non-vanishing. 

Severe limitations come also from baryogenesis considerations if one requires 
that the GUT scale baryon asymmetry is preserved. It has been shown in [15] 
that any bound from cosmology can be avoided by demanding that one of 
the lepton numbers is conserved over cosmological time scales. Since we wish 
to study the /x-e conversion, we could have conservation of r lepton number. 
On the other hand, if electroweak baryogenesis is assumed (see, e.g. [16]), the 



restrictions on lepton number are removed. Therefore we do not impose any 
extra constraints on lepton number violating couplings. 

Upper bounds on the magnitudes of individual couplings in Eq. (1) have been 
determined from a number of different sources. In particular one has consid- 
ered charged current universality, e — jj, — t universality, forward-backward 
asymmetry, z/^e scattering, atomic parity violation [17], neutrinoless double 
beta decay [18], v masses [13,19], Z-boson partial width [20], D^ , r-decays 
[21], D^ - Tf mixing, K+ , t-quark decays [22], hi production at LEP II [23] 
and heavy nuclei decay and n — n oscillation [24]. The bounds from these works 
have been collected and updated in recent reviews [25]. 

We should also note that very strong constraints on individual couplings A' 
can follow from the bounds on electric dipole moments of fermions [26] . How- 
ever, these bounds are qualitatively different from the previously mentioned 
ones, since they depend strongly on the existence of additional sources of 
CP-violation in the MSSM and the strength of the additional CP-violation. 
In fact, the upper bounds obtained from electric dipole moments apply for 
lA'Pl^l sin</), where 1^41 is the so-called A-term and is its phase, rather than 
for A'. If one makes the assumption that A is real and the only source of 
CP-violation in the MSSM is the CKM matrix then these limits are identi- 
cally washed away. Since these bounds are very model dependent and easily 
avoidable we shall not use them in our work. 

More recently there has been a lot of interest in deriving bounds on the prod- 
ucts of two such couplings which are stronger than the products of known 
bounds on the individual couplings. In addition to giving rise to unobserved 
rare processes these products turn out to be interesting also from the model 
building point of view. For example, bounds on these products constrain 
considerably models which attribute the origin of i?-parity violation to ap- 
proximate f/(l) [27] or f/(2) [28] flavour symmetry. Products of the type AA 
have been constrained from searches for muonium-antimuonium conversion 
[4], /i — i> 67 [29] and heavier lepton decays to three lighter ones [30]. Combi- 
nations of AA' have been bounded from the tree-level /i-e conversion, r- and 
TT-decays [4], from K ^ I1I2 [30] and from B ^ I1I2 [31]. Upper limits on A'A'-s 
are found from neutrinoless double beta decay [32], /x — > 67 [29][^ tree-level 
/i-e conversion [4], Am^, Attib [30], exotic top decays [33] and non-leptonic 
[34] and semi-leptonic [35] decays of heavy quark mesons. 

We have presented the combinations of AA and A'A' contributing to the /i-e 
conversion in Table 1 and Table 2, respectively. The upper limits on these 
combinations of couplings from the previous works are presented in the ta- 



^ The bounds on A'A'-s in Ref. [29] should all be divided by 3 due to the omission 
of the colour factor in the calculation. 
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fi-e at tree-level/ \/S 


/i-e 


at loop level/ ^/B 
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/ ~ 
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100 GeV 


100 GeV 


100 GeV 


1 TeV 


1 \' \' 1 

1^211 ^llll 


5.0-10-6 


8.0 


10-8 


1.1 


10-5 


8.3 - 10-4 


■^212-^112 


1.4-10-^ 


8.5 


10-8 


1.1 


10-5 


8.3 - 10-4 


1 \' \' 1 
1^213 ^113l 


1.4-10-^ 


8.5 


10-8 


1.1 


10-5 


8.3 - 10-4 


■^221 -^111 


5.0-10-6 


4.0 


10-7 


3.7 


10-5 


2.5-10-3 


1 \' \' 1 
1^222^^1121 


4.8 - 10-5 [30] 


4.3 


10-7 


3.7 


10-5 


2.5-10-3 


1 \' \' 1 
1^223^^1131 


4.8 - 10-5 [30] 


4.3 


10-7 


3.7 


10-5 


2.5 - 10-3 


K^231 -^111 


8.8-10-5 


1.0 


10-5 


1.2 


10-3 


8.3-10-2 


P^232-^112 


4.8-10-3 


1.1 


10-5 


1.2 


10-3 


8.3-10-2 


1 \' \' 1 
1^233 ^113 1 


4.8-10-3 


1.1 


10-5 


1.2 


10-3 


8.3-10-2 


■^211 -^121 


4.8 - 10-5 [30] 


4.0 


10-7 


7.3 


10-5 


2.5 - 10-3 


1 \' \' 1 
1^212^^1221 


4.8 - 10-5 [30] 


4.3 


10-7 


7.3 


10-5 


2.5 - 10-3 


1 \' \' 1 
1^213 ^123 1 


4.8 - 10-5 [30] 


4.3 


10-7 


7.3 


10-5 


2.5 - 10-3 


K^221 -^121 


1.4-10-^ 


8.0 


10-8 


2.0 


10-5 


1.2 - 10-3 


■^222-^122 


1.4-10-^ 


2.1 


10-6 


2.0 


10-5 


1.2 - 10-3 


1 \' \' 1 
1^223 ^123 1 


1.4-10-^ 


2.1 


10-6 


2.0 


10-5 


1.2-10-3 


^231 -^121 


2.6-10-3 


2.0 


10-6 


4.0 


10-4 


2.2-10-2 


P^232-^122 


4.8-10-3 


5.3 


10-5 


4.0 


10-4 


2.2- 10-2 


1 \' \' 1 
1^233 ^123 1 


4.8-10-3 


5.3 


10-5 


4.0 


10-4 


2.2-10-2 


K^211 -^131 


4.2 - 10"^ 
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10-5 


1.2 


10-3 


8.3-10-2 


■^212-^132 


4.0-10-3 


1.1 


10-5 


1.2 


10-3 


8.3-10-2 


■^213 -^133 


1.2-10-5 


1.1 


10-5 


1.2 


10-3 


8.3-10-2 


1 \' \' 1 
1^221 '^ISlI 


4.2 - 10-^ 


2.0 


10-6 
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10-4 


2.2-10-2 


1 \' \' 1 
1^222 ^132 1 


4.0-10-3 


5.3 


10-5 
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10-4 


2.2- 10-2 


1 \' \' 1 
1^223 ^133 1 


1.2-10-5 


5.3 


10-5 


4.0 


10-4 


2.2-10-2 


^231 -^131 


2.6 - 10-3 [29] 


8.0 


10-8 


8.7 


10-5 


3.2 - 10-3 


1 \' \' 1 
1^232 ^132 1 


2.6 - 10-3 [29] 


1.3 


10-3 


8.7 


10-5 


3.2 - 10-3 


1 \' \' 1 
1^233 ^133 1 


4.0-10-^ 


1.3 


10-3 


8.7 


10-5 


3.2 - 10-3 



Table 2 

Upper limits on the products jA'A'j for two different scalar masses m? =100 GeV 
and 1717=1 TeV. The previous bounds and the tree- level bounds scale quadratically 
with the sfermion mass. The scaling factor B is defined in the text and currently 
B = l. 5 



bles in the first columns. For A'A'-s the tree-level /i-e conversion bounds are 
presented separately in the second column of Table 2. All the couplings are 
assumed to be real. If the previously obtained bound on a particular prod- 
uct of couplings is stronger than the product of the bounds on the individual 
couplings then we have also presented in the tables the references from which 
the bound has been taken. In the opposite case the references are not given 
and can be found for example in Ref. [25]. The tree-level /x-e bounds are all 
taken from Ref. [4] assuming that there are no cancellations between different 
contributing terms. The bounds are given for two values, 100 GeV and 1 TeV, 
of the relevant supersymmetric scalar mass. If the listed bounds depend on the 
mass quadratically, which is the case for the bounds derived from tree-level 
processes, we have presented them for 100 GeV only. Extrapolation to 1 TeV 
can be done trivially by scaling the bound by two orders of magnitude. The 
quadratic dependence on the mass is also approximately true for the AA-type 
couplings in the loop-induced /i — > 67, but not for the A'A'-type couplings in- 
volving the third family. In this case the upper bound for rrif = 1 TeV becomes 
|AkA'l3.l<l•9■10-^ ^ = 1,2 [29]. 

The most stringent constraint on the products IA122A211I, |Ai32A3]^;^|, |Ai2iA']^;^2^| 
and |A23iA3^]^| follows from the tree-level fi-e conversion and is 4.0 ■ 10^^ for 
m^' = 100 GeV. 



3 Coherent fi—e conversion 



Due to the extraordinary sensitivity of ongoing experiments searching for /i-e 
conversion in nuclei, it can be expected that in the nearest future this process 
will become the main test of muon flavour conservation [5]. To take full ad- 
vantage of the /i-e experiments in studying possible new physics we present a 
general formalism for the coherent /i-e conversion. In practice only the appro- 
priate expressions for the coherent transition are needed, since the coherent 
process is enhanced by the large coherent nuclear charge and it dominates over 
the incoherent transition. Indeed, for the nuclei of interest the coherent rate 
of the conversion constitutes about 91-94% of the total rate [36]. 

The theory of //-e conversion in nuclei was first studied by Weinberg and 
Feinberg [37] and the early developments in this field have been summarised 
in [38]. Since then various nuclear models and approximations are used in the 
literature to calculate the coherent /i-e conversion nuclear form factors. It is 
important to notice that the results from the shell model [39], local density 
approximation [36] as well as quasi-particle RPA approximation [40] do not 
differ significantly from each other for both 22^^ and "^Ph nuclei showing 
consistency in understanding of nuclear physics involved [40]. We follow the 
notation of Ref. [36] taking into account the corrections to the approximation 



from the exact calculations performed in the same work. The corrections are 
negligible for 22^^ as the local density approximation works better for light 
nuclei but they are sizable, of the order of 40%, for "^Ph and Y^Au. 

The relevant /x-e conversion effective Lagrangian can be expressed as 

4vrQ; .;^ (pft) ^ Y^ / -A ji(n(m) 



_ ^^"-^ -A jKph) I "-^ V^ / -A ji(non) ■ jm{non)\ /r,\ 

'-eff — 2 J(vh)'^\ "•" fK Z^\Ji(non)-J\ -T Jm{non)-J ) , l^J 



i,m 



where the first term describes the photonic conversion and the second set of 
terms describe different contributions to the non-photonic conversion mecha- 
nisms. Here q^ denotes the photon momentum transfer, G the effective cou- 
pling constant appropriate for the interactions considered, j^ and j represent 
the leptonic and J^ and J hadronic vector and scalar currents, respectively. 
The photonic mechanism is enhanced by small q^ but can occur only at loop 
level and is thus suppressed by loop factors. The non-photonic mode is of inter- 
est, if the conversion can occur at tree-level like in models with non-diagonal 
Z' [41] and Higgs [42] couplings, models with leptoquarks [43], or if the loop 
contributions are enhanced by some other mechanism, e.g., models with non- 
decoupling of massive neutrinos [44] . In models with broken i?-parity, one may 
expect that both modes are important since the photonic mode is enhanced 
by large logarithms. 

Most generally the leptonic current for the photonic mechanism can be 
parametrised as 



.-A 
J(ph) 



Ueo + 75/Mo) 7. ig'-'' - ^\ + (/mi + 75/Ei) I cr'"^ 



/i(,3) 



where the form factors /eo, Iei, /mo and /mi can be computed from the 
underlying theory. Since the interaction of photon with fermions is well known, 
the hadronic current can in this case be written down immediately in terms 
of the nucleon spinor ^ = (p, n) as [38] 

^r = ^7A^(l + r3)^, (4) 



where ts is the isospin Pauli matrix. 

We do not consider non-photonic conversion mode in this Letter. However, 
since in the context of the /i-e conversion only the vector currents and the 
corresponding nuclear matrix elements have been discussed in the literature 
(see, e.g., Ref. [41]) while most of the models contain also new scalar particles, 
we feel that the following comment is in order. To calculate the tree-level /x-e 



conversion rate, one has to proceed from the quark level to the nucleon level by 
computing the quark current matrix elements (A^| q'yxq \N) = Gy N'jxN and 
(A^l qq \N) = G^' NN for the nucleon N =p, n (the coherent transition rate 
depends only on the quark vector or scalar couplings). Due to the conservation 
of the vector current and its coherent character, with the vector charge equal to 
the quark number, and assuming strong isospin symmetry, one can determine 
in the limit q^ ^ G^^'"^ = G^^'""^ = G^ = 2 and G^^'''^ = G^""^ = G^^^ = 1. 
However, in the case of the scalar operator qq one has to rely on nucleon 
models. Here we point out two extreme situations. In the limit of the fully non- 
relativistic quark model, for which qq can be basically regarded as a number 
operator, it is known that Gs ~ Gy- This also determines the upper bounds 
for the scalar form factors. The largest differences between the form factors 
are expected to occur in the case of fully relativistic models like the MIT 
bag model. In the latter case one has estimated the scalar form factors for 
massless quarks to be [45] Gg = 4/3 and Gg = 2/3. The actual values of 
the scalar form factors thus lay in between these two limiting cases. Therefore 
the constraints derived from the tree-level /i-e conversion in Ref. [4] can be 
somewhat relaxed due to this uncertainty but not more than a factor of 2/3 
or so. 

Due to the large mass of muon it is appropriate to take the customary non- 
relativistic limit for the motion of the muon in the muonic atom. In this limit 
the "large" component of the muon wave function factorizes out when cal- 
culating the coherent conversion rate. Therefore, for the photonic mechanism 
the coherent /x-e conversion branching ratio R^^ can be expressed as 



fl?]=c^p.£.9^z■^{„^ (5) 

q ^ capt 



where E^ is the electron energy. Eg ^ pe ^ m^ for this process, Tcapt is the 
total muon capture rate, |M(pe)P is the nuclear matrix element squared and 

eo = l/i?0 + /Ml|'+|/m + /Mo|' (6) 



shows the process dependence on the photonic form factors. The expression 
for |M(j9e)P in the local density approximation, 

mq)? = '^^^mq)\\ (7) 



the correction factors G to the approximation (compared with the exact cal- 
culation) as well as all numerical values of the above defined quantities for 
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22 



Ti and l2^Pb are presented in Ref. [36,46]. The result reads 



Rt = C 



8a'mlZ^ff\F,ip,)\'Z ^o' 



/^e ^ -p 4 

-*- capt y 



where C^' = 1.0, C^'' = 1.4, Zj/^ = 17.61, Z^J'f = 33.81, T^^^ = 2.59 ■ lO^^ 
s~^, r^pj = 1.3 ■ 10"^ s~^ and the proton nuclear form factors are Fp (q) = 

0.55 and Fp (g) = 0.25. Currently SINDRUM II experiment is running on 
gold, 79^ Am, but gold is not explicitly treated in Ref. [36]. However, since 
^Au _ ^Q gj^^ 2jPh _ g2 g^pg gQ close to each other then, within errors, all the 
needed quantities for ^^ Au and "^Ph are approximately equalQ. This result is 
strongly supported by theoretical calculations and experimental measurements 
of the total muon capture rate of Ph and Au [46] . In the numerical values we 
use the same experimental and theoretical input for both ^^ Au and "^^Ph. 



4 Bounds from the photonic /i— e conversion 



The /i-e conversion in supersymmetric theories which conserve i?-parity has 
been studied before in several works [47]. They conclude that other lepton 
flavour violating processes like /x — *■ 67 are about an order of magnitude more 
sensitive to the new phenomena than /x-e conversion. However, in the case of 
broken i?-parity one may expect the opposite result due to the logarithmic 
enhancement of /i-e conversion, which is about one order of magnitude at the 
amplitude level. 

Addition of Eq. (1) to the MSSM superpotential leads to more interactions in 
the model, but the MSSM particle content remains. The relevant part of the 
Lagrangian is found by standard techniques [48] 



+ '^ijk y^iLdkRdjL + djidkR^iL + dlfi^ii^djij 

-KJp (eidkRUpL + UpLdkRdiL + dkR^^iLUpL) + h.c, (9) 

where we have taken into account the flavour-mixing effects in the up-quark 
sector in terms of the CKM matrix elements Kjp. The effects of possible non- 
alignment of fermion and sfermion mass matrices are severely constrained [49] 
and therefore neglected. 



We thank H.C. Chiang and E. Oset for clarifying us this point. 
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Fig. 1. Contributing loop diagrams for yu-e conversion. Unless otherwise indicated 
i, A; = 1 ... 3. 



The Feynman diagrams contributing to the /i-e conversion at one-loop are 
shown in Fig. 1. The photon line is not shown, but it should be attached in 
all possible ways to the graphs. Following the general effective Lagrangian 
analyses of Ref. [5] we can say immediately that the first two diagrams with 
leptons in the loop give a logarithmically enhanced contribution to the /x-e 
conversion while the second two involving neutrino as an intermediate state 
are not enhanced. This is because in the latter set of diagrams one cannot 
attach the photon to the fermion line. For the typical fermion and sfermion 
masses we have |ln(mj/m|)| ~ (9(10) and, therefore, the contribution of 
the second set of diagrams can be neglected without affecting our numerical 
results. For the rest of the diagrams we compute the photonic form factors 
defined in Eq. (3). Starting from the Lagrangian (9) and using dimensional 
regularization we obtain for the form factors proportional to 7a 



J EO — ^J A 
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'MO 
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fir = -fLt = - J^ -V I In -V + Ffir) - ^ 1 , (12) 

and for the form factors proportional to (T^j, 






/mi - TIei - j^^, ^ , (13) 



/Ml -Ji^i - 24(47r)2 mj-' ^ ' 

fLf = f'ET' = 0. (15) 



Here the superscripts f = i, f = u and f = d denote contributions from 
the diagrams with leptons, u- and rf-quarks in the loop, respectively, r = 
m'j/{—q'^), where / stands for the virtual fermion in the loop, A'^c = 3 is the 
number of colours and 



Vl + 4r + r 



Ff{r) = — + 4r + Inr + (1 - 2r) Vl + 4r In \ . (16) 

3 \Vl + 4r - 1/ 

There are three important limiting cases. If r — ;► then Fy = — 1/3, if r ^ 1 
then Ff ^ 1.52 and if r > 1 then F/ = Inr + 4/3. 

As expected the form factors (10), (11), (12) contain large logarithms while 
the form factors (13), (14), (15) do not. The dependence of the form factors 
on particular combination of -R-parity violating couplings as well as on the 
± signs in the form factors coming from the leptonic loops should be fixed 
when calculating the constraints on the products of couplings. Deriving these 
expressions we have taken into account that g^ < and made a simplifying 
approximation Xf = m'j/m'^r <^ 1 which allows us to keep only the leading 
terms in Xf. This approximation is not adequate only in one case, i.e., if the 
fermion in the loop is top quark and 771? = 100 GeV. However, in this case 
the /i-e conversion is completely dominated by the diagram with the bottom 
quark in the loop and the top contribution can be neglected without making 
significant numerical error. Note also that in the first approximation Eq. (15) 
vanishes. For the expressions of /mi,_bi in higher orders of Xf we refer the 
reader to Ref. [29]. 

Now we are ready to constrain i?-parity violating couplings. Substituting the 
numerical values of the parameters to Eq. (8) we find 

RPjl = 8.6 (12.3) ■ 10"^ TeV^ ^ , (17) 
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where the first number in the expression corresponds to the conversion in Ti 
and the number in the brackets to the conversion in Pb and Au and ^q is 
defined in Eq. (6). Note that the conversion is somewhat enhanced in Pb and 
Au (in fact, maximised [36]) if compared with Ti. Substituting the current best 
bound on the conversion branching ratio Rj^l ^ 4.3 ■ 10^^^ [6] and the form 
factors above to Eq. (17) we find the constraints on the products of i?-parity 
violating couphngs from fi-e conversion at one-loop level. The bounds on the 
contributing combinations of AA are listed in Table 1 and on the combinations 
of A' A' in Table 2 for two values of sfermion masses tti ? = 100 GeV and 1 TeV. 

In future new experimental data on /i-e conversion will be available. To take 
into account the future improvements in the experimental sensitivity we define 
a scaling factor B as 



B 



RtKpre sent bound) I 1-0, iiNu = Ti, 



R^^{future bound) 2.8 , liNu = Au . 



To get new bounds one just has to divide the bounds in Table 1 and Table 2 
by the corresponding v5- 



5 Discussion and conclusions 



The photonic /i-e conversion which occurs at one-loop level probes the cou- 
pling products of the type AA and A'A'. For the combinations of AA comparison 
with the previously obtained bounds in Table 1 shows that in three cases out 
of six our new bounds are more stringent. This is because the conversion is 
enhanced by large logarithms ln(mj/?72|). If SINDRUM II experiment will 
reach the expected sensitivity and show negative results one has to scale the 
bounds down by factors of vB^ = 5 and vB^ = 12 which correspond to the 
expected hmits i?^^ ^ 5 ■ 10~^^ and Rj^^ ^ 3 ■ lO"^"', respectively. In this case 
the /i-e conversion bounds in Table 1 will be more stringent than the previous 
bounds for all testable combinations of AA. This is an important conclusion 
since no that significant improvements of the other experiments is expected. 
Note that in /i-e conversion AA-s can be probed only at loop level. 

The A'A'-type couplings give also rise to the tree-level /i-e conversion. The 
bounds derived from this process are very strong for most products of the 
couplings but not for all of them. The reason is that for some combinations of 
the couplings, especially if the third family squarks are involved, their contri- 
bution to the /i-e conversion can be strongly suppressed by small off-diagonal 
CKM matrix elements. In the worst case the suppression factor can be as 
large as A^, where \w ~ 0.2 is the Wolfenstein parameter. As a result the 
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tree-level fi-e conversion bounds are usually more stringent than the bounds 
from the other processes except for three cases. The bound on IA223A133I ob- 
tained by multiplying the bounds on the individual couplings is stronger than 
the fi-e conversion bound (note also that the previous bound on |A2i3A']^33| 
is only marginally weaker than the tree-level fi-e result). Importantly, as the 
result of our calculation, the bounds on IA232A132I and |A233A'^33| derived from 
the loop induced /i-e conversion are more stringent than the ones from the 
tree-level /x-e conversion. With the new data all the /i-e conversion bounds 
on |A'A'|-s will be more stringent than the others. Since the B factors for the 
loop induced and tree-level conversions are the same then the last two entries 
in Table 2 remain most stringently constrained by the loop results. 

In conclusion, searches for the /z-e conversion in nuclei, taking into account 
both the photonic conversion mode studied in this Letter and the non-photonic 
mode, constrain most of the products of -R-parity violating couplings testable 
in this process more stringently than do the other processes. In the nearest 
future this conclusion is expected to apply for all the products. 
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